Swelling of phospholipid floating bilayers: the effect of chain length 
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The equilibrium distance between two lipid bilayers stable in bulk water and in proximity of 
a substrate was investigated. Samples consisted of a homogeneous lipid bilayer, floating near an 
identical bilayer deposited on the hydrophilic surface of a silicon single crystal. Lipids were saturated 
di-acyl phosphocholines, with the number of carbon atoms per chain, n, varying from 16 to 20. The 
average and r.m.s. positions of the floating bilayer were determined by means of neutron specular 
reflectivity. Samples were prepared at room temperature (i.e. with the lipids in the gel phase) 
and measurements performed at various temperatures so that the whole region of transition from 
gel to fluid phase was explored. Data have been interpreted in terms of competition between the 
interbilayer potential and membrane fluctuations and used to estimate the bending rigidity of the 
bilayer. 

87.16.Dg: Membranes, Bilayers and Vesicles. 
87.15.Va: Fluctuations 
61. 12. Ha: Neutron reflectometry 
68.15.+e: Liquid thin films 

I. INTRODUCTION 

The basic building blocks of cell membranes are lipid bilayers. They provide mechanical stability and have a 
strong tendency to form closed structures. They support membrane proteins and have sufficient flexibility for vesicle 
budding and membrane fusion [1]. The physical properties of lipids are fundamental for understanding their biological 
functions. Because of their complexity the duplication of cell membranes and the investigation of their interactions 
with peptides or proteins is difficult. An increasing need for biomimetic models has recently caused the revival of a 
classical approach: the use of lipid bilayers [2] . A considerable effort has led to numerous successes in preparing new 
types of samples and corresponding structural measurement techniques. These can resolve sub-nanometric details, 
leading to investigations of lipid-lipid, lipid-peptide or lipid-protein interaction mechanisms. 

Phospholipids present a variety of different phases as a function of temperature [3]. Cooling from high to low 
temperature, lipids overcome the fluid to gel phase transition at a temperature commonly known as T TO (melting 
temperature). While in the fluid phase, L ai lipid chains are in a liquid- like conformation and mobile in the lateral 
direction, in the gel phase, Lp, chains are stiffer. For some lipids, just below T m , there is the so-called ripple 
phase, Pp, where the lamellae are deformed by a periodic and static out of plane modulation. Biophysical studies 
of membrane-membrane and membrane-protein interactions require well controlled model systems [2,4]. In previous 
papers [5,6] a new model system has been described consisting of stable and reproducible double bilayers, in which a 
lipid bilayer floats at 2-3 nm on top of an adsorbed one. Since the second bilayer is only weakly bound, the preparation 
is delicate but results are reproducible. The process relies on three vertical Langmuir-Blodgett depositions, followed 
by a horizontal Langmuir-Schaeffer one; this method is efficient when lipids arc in the gel phase and by raising the 
temperature the floating bilayer overcomes a phase transition becoming fluid and still stable. 

The floating bilayer system has been very useful to probe bilayer-bilayer interactions, and is currently being im- 
proved to accommodate for charges in water and in the bilayer, so that lipid/peptide interactions can be studied 
in physiological buffers. Since it has the potential to support transmembrane proteins, this system can be a better 
biomimetic model membrane than supported single bilayers. Floating bilayers offer the following advantages: they 
are well-defined structures, where the composition of each leaflet of the bilayer can be chosen separately; they are 
immersed in aqueous solution and they allow information on a single bilayer to be obtained [7,8], contrarily to the 
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classical models of multilamellar vesicles or stacked systems. They are stable in the fluid phase and swelling of the 
water layer between the bilayers can be observed with this system around the phase transition. 

A systematic study of the swelling is presented here, where the values of the distance of the floating bilayer from 
the adsorbed one were determined for phosphocholine lipids with saturated double chains and a number of carbon 
atoms, n, per chain from 16 to 20. For all samples a more or less broad region around the main phase transition 
temperature has been detected characterised by a swelling of the layer of water between the bilayers of about 1 nm. 
Giant swelling (i.e. swelling 3> lnm) was also observed with samples having 17 and 18 carbon chains. Results from 
a self-consistent theory [9] , that has been developed to extract the values of the bending modulus of the membranes 
from the experimental data, will be shown. 

II. MATERIALS AND METHODS 

Lipids were purchased from Avanti Polar Lipids (Lancaster, Alabama, USA) and used without further purification. 
They are phosphatidyl-cholines (PC) with different chain lengths (for simplicity they are indicated as Di-C„-PC): 
di-palmitoyl (Di-Ci6-PC), di-heptadecanoyl (D1-C17-PC), di-stearoyl (Di-Cis-PC), di-nonadecanoyl (Di-Cig-PC) and 
di-arachidoyl (Di-C2o-PC). Fully deuterated (Di-Cis-PC) (d83) was also used. Several attempts to prepare floating 
bilayers from fully deuterated Di-Ci6-PC failed. Hydrogenated lipids were measured in D 2 and fully deuterated lipids 
in H2O and Silicon Match Water, (i.e. mixture of H2O and D2O with a scattering length density = 2.07xlO~ 6 A -2 ). 
The minimum and maximum temperatures investigated are: 25.0-82.6°C for n=16; 25.0-55.6°C for n=17; 25. 0-64. 1C 
for n=18; 25.0-60.8°C for n=19 and 25.0-68.6°C for n=20. 

Two kinds of silicon substrates were used: (i) of surface 5x5 cm 2 and 1 cm thick obtained and polished by SESO 
(Aix-en Provence, F) with a r.m.s. roughness 0.15 nm measured on 300x300 /1m 2 ; and (ii) of surface 8x5 cm 2 and 2 
cm thick obtained from Siltronix (Archamps, F) and polished by the Optics Laboratory of the European Synchrotron 
Radiation Facility (Grenoble, F) with a technique including diamond powder treatment and producing high quality 
crystal surfaces with r.m.s. roughness of <0.2 nm. 

The silicon blocks were made highly hydrophilic with a UV/ozone treatment as in [6]. Lipids were deposited on 
the substrate with a combination of the Langmuir-Blodgett and Langmuir-Schaeffer techniques as described in [5,6]. 
Depositions were done at room temperature, and the lipids compressed to surface pressures varying with the lipid 
chain length (40 mN/m for n=16,17,18; 35 mN/m for n=19 and 30 mN/m for n~20). Transfer ratios for the first 
layers were always high and of the same order as samples described in [6] . They improved as the chain length increased 
and attained the value of 1.0 for n=18,19,20. It was found that the horizontality of the sample in the Schacffcr dip is 
crucial for a good quality sample. In fact, while in [5,6] during the Schacffcr dip the substrates were slightly inclined 
so that no air bubbles would form on the solid surface, here systematic measurements on Di-Cis-PC samples were 
made and it was found that a sample surface as parallel as possible to the water surface (micrometric screws were 
used) prevents the formation of bubbles and gives transfer ratios always bigger than 0.97. 

While it is difficult to prepare samples with short phospholipid chains (e.g. with DMPC), it has been found that 
this is possible if the first layer is chemically grafted to the surface [10], although in this case the coverage of the 
DMPC floating layer is not always optimal and reproducibility problems are encountered. This method when applied 
to longer chain lipids has resulted in the formation of very stable and defect-free layers (work in progress). 

III. NEUTRON REFLECTIVITY: BASIC PRINCIPLES AND MEASUREMENTS 

Specular reflectivity, R(Q), defined as the ratio between the specularly reflected and incoming intensities of a 
neutron beam, is measured as a function of the wave vector transfer, Q = AirsmO/X perpendicular to the reflecting 
surface, where 8 is the angle and A the wavelength of the incoming beam. R(Q) is related to the scattering length 
density across the interface, p(z), and for values of Q bigger than the critical value for total reflection, it decays as 

Q- 4 [11]. 

Specular reflectivity allows the determination of the structure of matter perpendicular to a surface or an interface. 
Experiments are performed in reflection at grazing incidence. Samples must be planar, very flat and with roughness 
as small as possible (<0.5 nm). The data presented here were collected at the High Flux Reactor (HFR) of the 
Institut-Laue Langevin (ILL, Grenoble, F) both on the small angle diffractometer D16 and the high flux reflectometer 
D17. Instrument set-up and the data collection method on D16 have been described in [5]. A monochromatic neutron 
beam of wavelength A= 4.52 A impinged on the sample at grazing incidence. A two-dimensional 3 He wire detector 
was kept in a fixed position covering the angular 29 range 0°- 9°. The beam was vertically focused at the detector 
position in order to gain intensity without worsening divergence. The dual mode instrument D17 [12,13] is designed 
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to take advantage of both Time-Of-Flight (TOF) and monochromatic methods of measuring reflectivity. In this study 
measurements were all taken with the instrument in the TOF mode by using a spread of wavelengths between 2 and 
20A at two incoming angles (typically 0.7° and 4°) and with a time resolution, defined by two rotating choppers, of 
At/t = 1% and 2-5%, respectively, where At is the neutron pulse duration. The beam was defined in the horizontal 
direction by a set of two slits, one just before the sample and one before a vertically focusing guide. Neutron reflectivity 
from the same sample was measured at different temperatures monitored with a thermocouple (equilibration time 
25 minutes, stability < 0.1°C, absolute precision < 0.3°C), in the water-regulated sample chamber already described 
in [5]. Two types of measurements were performed: (i) up to Q ~ 0.08 A -1 ; (ii) over the whole available Q-range 
defined as follows. 

For measurements done on D16, the whole available Q-range (0.003-0.22 A" 1 ) was scanned in eight hours, after 
equilibration. On D17 the useful Q-range, before hitting the sample background, spanned from 0.007 to 0.25 A -1 
and was measured in about two hours, also after equilibration. The equilibration time varied depending on the step 
of temperature increase or decrease. For a step of 1°C it was of about 15 minutes. 

IV. DATA ANALYSIS: MODEL FITTING OF REFLECTIVITY DATA 

A reflectivity curve essentially reflects the sample density perpendicular to the substrate surface, or rather the 
square modulus of its Fourier transform. Data analysis requires a good knowledge of the sample and here it was done 
via model fitting. Born and Wolf give a general solution, the so called optical matrix method [14], to calculate the 
reflectivity from any number of parallel, homogeneous layers, which is particularly useful since any layer structure 
can be approximately described by dividing it into an arbitrary number of layers parallel to the interface, each having 
a uniform scattering length density. In previous studies [6,15], multiple contrast neutron measurements at room 
temperature had allowed the determination within A precision of the profile of adsorbed and floating bilayers. Each 
bilayer is resolved into outer-inner-outer slabs: outer slab = heads (phospho-choline and glycerol groups), inner slab 
= tails (hydrocarbon chains). A water film of thickness T) w separates the adsorbed and floating bilayers. Together 
with the natural oxide layer present on the silicon crystal and the hydration water layer between the substrate and 
the first lipid layer, this leads to 9 slabs in total, each one having its own average thickness and scattering length 
density. Each interface between two slabs has a certain width, also called r.m.s. roughness. Since this is included into 
the density profile, it can be determined by fits of reflectivity curves, but it is not possible to determine whether it is 
a static intrinsic width or an average over the neutron beam correlation length of temporal or spatial fluctuations. 

The present work focuses on thicknesses and roughnesses only. By starting from the above 9-slab model, with the 
tabulated values of p (see Tables in refs. [5,6]) the relevant parameters to the data were fitted for all systems in the gel 
phase (the starting point of all measurements was the collection of the reflectivity profile at 25°C). The first reflectivity 
minimum abscissa is very sensitive to the thickness of both water films. The second reflectivity minimum position is 
mainly sensitive to the thickness of the first bilayer tail slab; the amplitudes and sharpnesses of both minima depend 
upon bilayer roughness (sec figures 1 and 2). All slabs were easily distinguished, except that the high hydration of 
lipid heads [16,17] prevents a localization of the heads/D20 interface to better than 0.1 nm. Correcting the scattering 
length densities for temperature variations had no significant effect on the positions of the reflectivity minima, hence 
on the fitted parameters. For data at high temperatures the analysis concentrated on the exact determination of T) w . 
When giant swelling occurred, data analysis was complicated by the fact that the high values found for the roughness 
limit the usefulness of the analysis based on homogeneous layer models, leading to large error bars for D w and a. 

In a parallel study (in progress) of specular and off-specular reflectivity from synchrotron radiation on similar 
samples the wider q-range attained and therefore the inadequacy of the box model fitting [18] has made it necessary 
to use the one Gaussian hybrid model described in [19]. The results as far as D w is concerned are in perfect agreement 
with the values obtained from neutron reflectivity using the box model. 

Figure 1 gives some examples of the collected data. The reflectivity curves obtained from hydrogenated lipids in 
the fluid phase at randomly chosen temperatures, as well as the best fits to the data, are shown. As expected the 
main differences are observed at the second minimum which is dependent on the single bilayer structure. Figure 2 is 
an example of curves collected for the system Di-Ci7-PC in the gel and fluid phases and in the transition region, as 
well as the best fits to the data. The density profiles obtained from the model that best fits the data are given in the 
insert. 
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V. DATA ANALYSIS: CALCULATION OF BENDING MODULUS 



The self-consistent theory described in [9] for a fluid bilayer floating near a solid substrate was applied to the results 
of this study. The basic principles are briefly described below. The shape of the substrate potential determines the 
amplitude a of thermal fluctuations which itself determines the effective potential [3,20] felt by the bilayer. The 
average bilayer position and its r.m.s. fluctuation amplitude are thus coupled and depend both on temperature and 
bending energy of the membrane. This model situation is experimentally relevant as already reported in [4]. Three 
self-consistent equations may be written for the partition function, Z, the equilibrium inter-bilaycr distance, D^, and 
the r.m.s. fluctuation amplitude, a: 

• , -16 „ K s , er 2 U(z)-3 (z ~ D ^ )2 

Z = dze < fc B T >" 8 " 2 , 



D w = IJdzz e - 16 TW^- 31 ^, (1) 
a* = Udz(z- D w f e - 16 cW^>- 3 ^ . 



where fee is the Boltzmann constant, T the temperature and k the bilayer bending rigidity modulus. U(z) is a one 
dimensional potential describing the interaction between the bilayer and the substrate. 

At short distances, with z typically of the order of the hydration length zq ~ 0.6 nm or smaller, the potential is 
dominated by the hydration repulsion due to water molecules inserted between hydrophilic lipid heads [3] : 

U h (z) = A h exp (-2^) • (2) 

At long distances, for a neutral membrane, the dominant term in U(z) is due to van der Waals attraction : 

f 2 1 



U v dw{z) = -7^— 

1Z7T 



z 2 {z + Sf (z + 25) 2 



(3) 



where S is the bilayer thickness and A the Hamaker constant. 

Here, it is assumed that A does not diverge and does not undergo any marked minimum so that its variations with 
temperature are negligible (compared to those of k) and therefore A is constant in what follows. 

The dimensionless parameter (3 : 

determines the position and the fluctuations of the membrane. At low values of (3, the bilayer is stiff, barely 
fluctuating in a strong minimum of potential energy. At high values of f3, the bilayer is soft and has large fluctuations 
leading to the increase of the mean distance D w to the adsorbed bilayer. 

By eliminating Z from the three self-consistent equations (1), two of the three physical quantities D, a and (3 can 
be expressed as a function of the third one. Experimentally, D„, and a are measurable quantities while j3 is much 
more difficult to measure, especially for bilayers in the gel phase [21]. Therefore, measurements of D w are used here 
to estimate (3 and a by using a classical numerical resolution of equations (1) [9]. The validity of this approach is 
tested by comparing a to the experimental data. 



VI. RESULTS AND DISCUSSION 



The analysis of the reflectivity curves has led to the determination of the water layer D^, and the bilayer roughness 
<7 (see section IV). The main results are summarised in Table I as well as in Figures 3 to 9. It was found that the 
silicon crystals were covered by a 1.2±0.1 nm thick oxide. A thin water film of thickness 0.6±0.1 nm separated it from 
the adsorbed bilayer. Both in the gel and in the fluid phases the thickness of the headgroups was 0.8±0.1 nm. The 
thickness of chains (D c hain) in the gel phase (at 25°C) is reported in Table I for all lipid species. In the fluid phase 
those thicknesses were usually lower by about 0.2 nm. The bilayers were found to incorporate 5 to 10% of water. As 
for the roughness, well below and above the transition region it was the same as the substrate, i.e. 0.3±0.1 nm. In 
the transition region it increased (see below). 
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Figures 3 and 4 are a summary of the changes of D w with temperature for all lipid species. Figure 3 refers to 
measurements from hydrogenated lipids in D2O and Figure 4 refers to measurements of deuterated Di-Cis-PC in 
H 2 and contrast match water. On the x-axis is T-T m , i.e the difference between the temperature at which the 
measurement was taken and the value of the melting temperature (the values used for T m are listed in Table I). For 
all lipids a very similar behaviour was observed, that is an increase of both D w and a in a certain range of temperature 
around the main phase transition (swelling region) . The degree of swelling was always around 1 nm, and in two cases 
(n=17,18) giant swelling (^>1 nm) was observed. The swelling temperature T s defined as the temperature at which 
D w and a are maximun appears to be systematically 4 to 5 degrees below the melting temperature (sec also figure 
9). Above the transition, the floating bilayer is stable and D w goes back to values similar to those found in the gel 
phase, although slightly higher. This equilibrium distance in the fluid phase increases slowly with temperature. In 
the gel phase, the minimum value D w _ g , listed in Table I, is obtained at the lowest measured temperature (i.e. 25°C). 
This value seems to be constant over a large temperature range, indicating that it is probably a good estimation of 
the bilayer position in the absence of thermal fluctuations. Figure 5 shows the variation of this minimum value with 
the lipid chain length n. At 25°C, D w , g tends to decrease when n increases, and the effect seems to depend also on 
the oddness or evenness of n. It is unclear if this is fortuitous, to the best of the authors knowledge there are no other 
examples in the literature. 

TABLE I. List of experimental data and literature parameters for the lipid species studied, n is the number of carbon 
atoms per chain in the formula Di-C„-PC; T p is the temperature of transition from the gel to the ripple phase [22], T m the 
temperature of main transition [22]. In column 4 are listed the regions of temperature where swelling was observed, T s is 
temperature at which D^, reaches its maximum value D m (max) and D w ,g is the value of in Gel phase. T> c hain is the chain 
region thickness, for the floating bilayer, in the gel phase (25°C). Ah and A are defined in the text and estimated from D„, s 
(see below). 



n 


T P 

(°C) 


T m 
(°C) 


Swelling Region 

(°C) 


T s (max) 
(°C) 


Deft, ain 

(nm) 


D„ (max) 
(nm) 


D-u;, g 

(nm) 


A 


16 


35-37 


41-41.2 
(41) 2 


(35-40) ±1 




3.2±0.2 


3.3±0.2 


2.2 ± 0.1 


80 ± 15 


17 


43 


48-49.8 
(48) 2 


(40-45) ±1 


41.5±0.5 


3.4±0.2 


4.8±0.2 


1.4 ± 0.1 


23 ± 3 


18 


50-52 


54.6-55 
(55) 2 


(45-55)±l 


51.5±0.5 


3.6±0.2 


6.0±0.4 


1.8 ± 0.1 


41 ± 7 


18 (D) 1 




(50.5) 


(37-50) ±1 




3.6±0.2 


2.8±0.1 


1.8 ± 0.1 


41 ± 7 


19 


55.7-57.8 


59-60 
(60) 2 


(50-62)±l 




3.8±0.2 


2.3±0.1 


1.2 ± 0.1 


19 ± 2 


20 


62.1-63.7 


63-66.4 
(66) 2 


(56-64) ±1 




4.0±0.2 


2.5±0.1 


1.5 ± 0.1 


26 ± 3 



1 Fully deuterated lipid; 2 Value given by Avanti polar Lipids and used for the figures; 



5 



The reproducibility of the observed effects was checked on different samples, from different experiments carried out 
in a period of two years on various silicon substrates and instruments (either D16 or D17). Figures 6(a) and 6(b) show 
results from samples with n=16 and 18. For clarity results have been split in different graphs so that the most notable 
effects could be highlighted. For the n=18 case data from 4 samples are reported (samples A, B, C, and D in the 
caption). The temperature range where the swelling was observed is reproducible. The amplitude of swelling, that is 
the maximal distance reached, is variable from one sample to the other. In two cases giant swelling was observed, i.e. 
with lipids having n=17 and n=18, the last one already described in [6]. It is not excluded that giant swelling occurs 
also with the other lipids investigated in this geometry The limited available beam-time has not allowed exploring 
the transition temperature ranges in fine detail for all samples and since giant swelling is usually observed in narrow 
T ranges [23,24] it might have been missed in other cases. 

The reversibility of the swelling was also studied. Figures 7(a) and 7(b) show results from one sample with n=17 
and two samples with n=18 with data collected while heating up and cooling down once or twice, as indicated in the 
figure caption. The major result is that the swelling is reversible. On both figures, a weak increase in the maximum 
amplitude of the swelling with the successive heating and cooling processes is observed. It could still be due to the 
narrowness of the temperature range where swelling occurs, as mentioned above or it could be an indication that some 
defects of the bilayer like pinning points have been annealed. The swelling temperature is barely sensitive to cycles 
in temperature : T s appears to be slightly lower when cooling down than when heating up. This could be related to 
the hysteresis of the main transition although the effect is too weak for a definitive interpretation. 

Ref. [9] suggests a possible interpretation of the swelling in term of a balance between attraction by the substrate 
and entropic repulsion due to thermal fluctuation. The self-consistent theory briefly described in section V enables 
expressing D m and a as function of (3 by assuming that the dimensionless parameter ^ and the characteristic length 
zq of the microscopic potential (equations 2 and 3) are known. Here it was assumed that zq is equal to 0.6 nm and 
that it does not depend on chain length. By using the experimental value of D ffi at low temperature (D w , g in table I) 
and the microscopic potential, the ratio ^ was estimated for different values of n (Table I). From this approach, the 
experimental values of T) w are used to estimate (3 and the corresponding value of a for each temperature. The value 
of a is then compared to the r.m.s. roughness obtained by model fitting the data. Results from one representative 
example (i.e. n=16) are given in figure 8. Agreement with all samples is remarkable (data not shown). 

The derivation of an experimental estimation of n implies the knowledge of the energy scale A (Hamaker constant). 
This value was estimated for all samples by normalising k to the usual value lOfc^T in the fluid phase. Values of A 
of the order of 10~ 21 — 1CP 22 J were found in agreement with literature data [3]. This method forces the values for all 
different lipids to overlap artificially at high temperature. Figure 9 shows the value of the bending modulus calculated 
in that way. A minimum of the bending modulus is found around the main phase transition. Such minimum in k is in 
agreement with literature data [25-28], even if its physical origin is unclear. It could be a signature of the proximity 
to a critical point in the phase diagram. As already mentioned, for Di-Cis-PC in D2O the pre-transition temperature 
T p = 51-52°C [29,30] is close to the temperature at which the largest swelling is observed here. On the other hand, 
a minimum in k is more likely to exist near the main (melting) transition T m [27,31]. This would imply that for the 
floating bilayer T m is lower than the value reported for multilamellar systems [29,30]. A comparable lowering of the 
melting temperature has been previously observed on single bilayers on a solid support [32] and was attributed to 
surface tension effects. Two experimental facts enforce this hypothesis: (i) the swelling temperature T s seems to be 
always shifted by 4-5°C from the bulk T m value (figure 3 and 9) independently of the phospholipid chain length; (ii) 
measurement made on deuterated Di-Cis-PC in H 2 in water exhibit also swelling and no literature data were found 
on a transition to the ripple phase for those deuterated species. In floating bilayers, a lateral stress might arise from 
possible pinning on defects. However, such a large tension would strongly limit the fluctuations and prevent both 
and a from increasing as much as we observe [33,34]. Finally, this interpretation can explain why the equilibrium 
distance is highly sensitive to small variations in the bending modulus and in the balance between attraction and 
repulsion (the samples are on the verge of unbinding) [33,35]. 

VII. CONCLUSIONS AND PERSPECTIVES 

Floating bilayers offer interesting advantages over other models for membranes. A structured bilayer, stable in time 
and free to fluctuate in the fluid phase enables in-situ studies of transmembrane channels or interactions between a 
lipid bilayer and membrane proteins. Its size, orientation, and the proximity of a smooth substrate enable reflectivity 
studies, and therefore structural characterisation in the fraction of nanometre sea le. The composition of both lipids 
and surrounding medium can be varied with a certain freedom. Furthermore, this is a well defined geometry for 
comparison with theoretical models; it also provides different control parameters (temperature, composition of lipids 
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or water) to study the interaction between two membranes. Thanks to the excess water, a large increase in both 
the floating bilaycr roughness and the inter-bilayer distance has been observed around T m for a series of samples 
from phosphocholine lipids having different chain lengths. Anomalous swelling has been observed frequently at the 
main phase transition in multilamellar systems and it is object of debate [36,37] whether it is due to an increase of 
the bilayer thickness or to the swelling of the water layer between the membranes. The results presented here can 
help in this debate since the swelling that we observe is undoubtedly due to the increase in thickness of the water 
layer between the two bilayers. A parallel study is in progress for the determination of the fluctuation spectrum of 
the floating bilayers by using synchrotron radiation off-specular reflectivity measurements. In fact, while specular 
reflection gives information in the direction perpendicular to the interface, the lateral structure of the interface may 
be probed by the nonspecular scattering measured at reflection angles different from the specular one. 

Since in that case the intensity of the reflected beam is much weaker than in the specular direction, the use of 
synchrotron radiation is necessary for obtaining a wide enough dynamic range for the analysis of these thin systems. 
Measurements have been performed both on double bilayers and on samples prepared by depositing three lipid layers 
on hydrophobically grafted silicon substrates. In this way the fluctuations of the layer closer to the substrate are 
reduced and information can be obtained only from the floating one. Results (to be published) indicate that the layer 
is under some tension (of the order of 3-5 mN/m) and such tension could explain why we observe the swelling at 
values of temperature lower than literature values of T m . 

Besides the bilayer fluctuation studies, the model is being improved to make it more relevant for interactions with 
biological systems. So far, we have succeeded inserting various amounts of cholesterol ranging from 0.5 to 10% choles- 
terol in the bilayers as well as negative charges (work in progress). The preparation of asymmetric bilayers, with the 
internal leaflet formed by phospho-ethanolamine headgroups and the outer one by phosphatidyl-choline headgroups 
is being optimised. 
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VIII. FIGURE CAPTIONS 



FIG. 1. Reflectivity data (points) and best fits to the data (continuous lines) from all lipid species, Di-C„-PC, in the 
hydrogenated form and measured in D 2 0, in the fluid phase: (•) n=16 at T=44.2°C, (o) n=17 at T=45.8°C, (■) n=18 at 
T=55.4°G (□) n=19 at T=58.9°C, (0) n=20 at T=68.6°C. 



FIG. 2. Reflectivity profiles (points) and best fits to the data (continuous lines) profiles for a double bilayer of Di-Ci7-PC 
at (•) T=25.0°C, (o) T=43.0°C and (□) T=45.8°G In the insert are the scattering length density profiles as from the model 
used to fit the data: continuous thin line for T=25.0°C, continuous thick line for T=43.0°C and dashed line for T=45.8°C 



FIG 3. Distance, D TO , between the two bilayers from measurements from hydrogenated lipids vs. the deviation of temperature 
from the melting temperature, T-T m ; (•) n=16; (o) n=17; (■) n=18; (□) n=19; (0) n=20. Lines between points are added 
for clarity. Samples with n=16 and 18 were measured on D16; samples with n=17,19 and 20 were measured on D17. 

FIG 4. Distance, D w , between the two bilayers from measurements from deuteriated Di-Cis-PC in Silicon Match Water (o) 
and H2O (A) and for hydrogenated Di-Gg-PC in D2O (■) (same datas as figure 3) vs. the deviation of temperature from the 
melting temperature, T-T m . 



FIG 5. Variation of D„ as a function of n at T=25°C. 



FIG 6. (a) Distance, D„, between two Di-Ci6-PC bilayers as a function of temperature, T. Experiments performed at 
different times on surfaces of different size: (o) 5x5 cm 2 (•) 5x8cm 2 . Both experiments were performed on D16.(b) Distance, 
D„, between two Di-Cig-PC bilayers. Experiments performed at different times over two years on different instruments: (□) 
sample A, measured on D17; (v) sample A, second time raised in temperature after cooling from 58°C; (o) sample B, measured 
on D16; (A) sample C measured on D16; (O) sample D measured on D16 [6]. 

FIG 7. (a) Distance, D„, between two Di-Ci7-PC bilayers as a function of the temperature, T: (□) first time up in T; (•) 
first time down in T; (O) second time up in T. (b) Distance, Dw, between two Di-Gs-PC bilayers for 2 samples. Sample 1: 
(□) first time up (1); (■) first time down (2); (O) second time up (3); Sample 2: (o) first time up (1'); (•) first time down (2'). 
All experiments were performed on D17. 

FIG 8. Comparison of values of the r.m.s. roughness for a Di-Ci6-PC sample heated from 25 to 82.5° C as a function of 
temperature, T (□) as determined experimentally by model fitting the data and (•) as extracted with the self-consistent theory 
[9]. 

FIG 9. Values of the bending modulus k normalised by fcsT as a function of the deviation of temperature from the melting 
temperature, T-T m for the samples Di-C„-PC with (•) n=16; (o) n=17; (■) n=18; (□) n=19; (<>) n=20. Values of k are 
normalised to the arbitrary value lOfesT in the fluid phase. This method forces the values for all different lipids to overlap 
artificially in the fluid phase. The dashed line represents the pretransition temperature T p . 
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